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I.  INTRODUCTION 


The  interior  ballistic  process  in  the  regenerative  liquid  propellant  gun 

is  primarily  controlled  by  the  rate  of  injection  of  the  liquid  propellant, 

and,  thus  by  the  (notion  of  the  regenerative  piston.  In  the  interior  ballistic 

1  2  3  4  5 

models  developed  to  date  by  Gough,  Coffee,  Cushman  and  Bulman  the 

equations  of  motion  for  the  regenerative  piston  describe  only  the  pressure  and 
friction  forces.  These  models,  in  general,  neglect  any  direct  coupling 
between  piston  motion  a..d  liquid  injection.  The  formulation  of  Cushman,4 
which  has  been  discussed  by  Gough, ^  is  an  exception. 

Using  Gough's^  approach,  the  equations  of  motion  for  the  piston  and  the 
liquid  can  be  written  in  the  form. 


The  first  two  terms  on  the  right  hand  side  of  the  equation  account  for 
the  kinetic  energy  of  the  liquid  exiting  and  approaching  the  injector, 
respectively.  The  term  involving  ip  is  an  injector  entrance  loss,  as  described 
by  Kaufmann,  and  the  last  term  is  the  usual  Blasius  correlation  for  the 
pressure  head  loss  due  to  friction  in  turbulent  pipe  flow. 


1 


As  noted,  the  more  usual  formulations  of  the  equations  for  piston  motion 
and  liquid  injection  neglect  coupling  between  the  physical  processes.  In 
general/  the  acceleration  of  the  liquid  through  the  injector  is  also 
neglected/  resulting  in  equations  of  the  form/ 

•  ’*■ 

U  =  1/M«  {PA  “  }  '  (4) 

p  p  3  p  R 


v3  =  CD  [2  <P  -  P3)/pL] 1/2  .  (5) 

While  it  is  acknowledged  that  the  development  of  the  above  equations 

involve  numerous  approximations/  until  recently  it  has  appeared  that  such 

simple  models  were  adequate  to  describe  piston  motion  and  liquid  injection  in 

the  regenerative  gun.  However,  Pate7  has  reported  values  of  the  discharge 

coefficient  in  cold  flow  experiments  in  a  regenerative  LP  fixture  which 

significantly  exceed  those  predicted  using  equation  3.  In  addition.  Coffee8 

has  reported  a  discharge  coefficient  derived  from  regenerative  gun  firing  data 

which  not  only  exceeds  the  predicted  value,  but  also  is  highly  transient  in 

g 

nature.  Finally,  Rizk  and  Edelman  have  reported  values  of  the  discharge 
coefficient,  obtained  from  a  two  dimensional  simulation  of  the  cold  flow 
experiment  of  Pate,  which  are  in  agreement  with  experimental  data. 

The  objective  of  this  work  is  to  investigate  liquid  injection  in  the 
regenerative  liquid  propellant  gun,  within  the  context  of  a  lumped  parameter 
model.  This  requires  the  development  of  a  set  of  ordinary  differential 
equations,  similar  to  1  and  2,  but  which  more  accurately  account  for?  (1)  the 
coupling  between  piston  motion  and  liquid  injection  and  (2)  the  inertia  of  the 
liquid  in  the  reservoir.  This  report  describes  the  progress  to  date  in  this 
effort.  A  general  mathematical  model  describing  the  process  has  been 
developed.  Results  of  the  simulation  are  presented  and  compared  to 

experimental  data  from  a  30mm  RLPG  test  fixture. 


The  governing  equations  have  been  coded  in  FORTRAN  for  the  IBM  PC-AT, 
using  the  Adams  method  with  functional  iteration  from  the  IMSL  Library,  DGEAR. 

II.  EQUATIONS  OF  MOTION 

The  equations  of  motion  are  written  for  the  control  volume  containing  the 
regenerative  piston  and  the  liquid  propellant  reservoir,  as  shown  in  Fig.  1. 


CONTROL  VOLUME 


o  CONTROL  VOLUME  INCLUDES 
RESERVOIR  AND  PISTON 


Figure  1.  Control  Volume 


•Hie  contours  of  the  piston  and  the  reservoir  are  approximated  by  straight 
line  segments  as  indicated.  The  center  bolt  and  transducer  bolt  are  fixed  in 
the  reference  frame  of  the  chamber.  The  origin  of  our  coordinate  system  is 
fixed  at  the  rear  (left  hand)  end  of  the  reservoir,  and  x  is  the  coordinate 
along  the  bolt  with  Xj»  x2»  and  x^  the  coordinates  of  fixed  positions  on  the 
bolt  as  shown.  The  piston  moves  rearward  with  a  velocity  u^,  and  the  points 
s^'  a2'  an<*  ®3  are  coordinates  of  fixed  stations  on  the  inner  contour  of 
the  piston  with  respect  to  the  origin,  as  shown,  such  that  these  coordinates 
vary  with  time  as  the  piston  is  displaced  to  the  left.  Note  that  the  right 
hand  face  of  the  control  volume  is  attached  to  the  chamber  face  of  the  piston 
( s 3)  such  that  the  control  volume  also  varies  with  time. 

The  equations  of  motion  are  written  in  the  reference  frame  of  the  bolt 
(chamber).  The  momentum  equation  for  the  control  volume  [0,83]  is 


where  Mp  is  the  mass  of  the  piston  and  dA  is  the  outward  directed  normal  from 
the  element  of  control  surface.  Then, 


~M  u  it 
p  p  3t 


I: 


p  vA  dx 

li 


i  +  PL  V3  A3  1 


[P0  \  +  PCF  AS  "  P3  (Ap  +  A3)]  1 


(7) 


where  the  control  volume  extends  to  include  the  piston  shaft  and  P 

L.r 

represents  the  pressure  of  a  control  fluid  on  the  area  of  the  piston  shaft, 
Ag.  We  will  ignore  this  term  in  this  paper,  and  write  the  momentum  equation 
for  the  control  volume  as 
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•*» 


n*  i 


$' 

I1 


The  unsteady  Bernoulli  equation  is 


(9) 


where  hf  is  the  headloss  due  to  friction,  and  h^’  is  an  entrance  loss. 
The  integrals  of  interest  are  then. 


v(x,t)  dx 


and 


PL(x,t)  v(x,t)  A(x,t )  dx 


In  order  to  evaluate  these  integrals,  we  require  p(x,t)  and  v(x,t)  . 


III.  LAGRANGE  APPROXIMATION  WITH  AREA  CHANGE 


The  equations  of  motion  for  the  fluid  are 


where  x  indicates  the  position  on  the  bolt,  R1  indicates  the  radius  of  the 
piston  at  s^  and  is  the  radius  of  the  piston  at  s^.  H(x)  is  the  Heaviside 
function  defined  as  H(x)  =  0  for  x<0  and  H(x)  =  1  for  x>=0.  The  radius  of  the 
inner  surface  of  the  piston  at  a  position  x  is  time  dependent  since  the  piston 
is  moving. 

Similarly,  the  radius  of  the  bolt  at  any  position  x  on  the  bolt  can  be 
expressed 
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rKrirt.rilSl*.  iVnIVwi 


rb(x)  *  {rb  +  -- — Z~x~)  [x  "  X1H1  -  H(x1  -  x)]}  H(x2 


+  r2  [1  -  H{x2  -  x)]  H(x3  -  x)  • 


The  cross-sectional  area  of  the  liquid  is  then  given  by 


2  2 

A(x,t)  -  n  [R  ( X/ t )  -  r  (x)] 

b 


and  the  volume  of  the  liquid  is 


V(x,t) 


A(x',t)  dx'  . 


Then  it  can  be  shown  that 


A(x,t) 


3A(x,t) 
p  3x 


and 


V 


R 


u 

P 


(A  +  A  )  . 
R  3 


Returning  now  to  equation  (10),  and  assuming 


3p  /3x  =  0 
L 


7 


We  have 


pl  at 


_3A_  _  _3vA_ 
3t  3x 


,PL  3  K1 

3t  '  5t  LvJ 


m  mV 

L  L  R 


-A- 

R 


PLV3A3  '  PLUpA3]  ”  ["PLUp  (AR  +  A3)] 


where 


pv^A^  is  the  mass  flux  through  the  orifice  with  respect  to  the  bolt,  and 


pu  A  is  the  additional  mass  flux  into  the  chamber  due  to  piston  motion. 
P  3 


Then, 


PL  3t 


V3A3  "  Vk 
. “ 


Now  using  (18)  and  (21),  equation  (20)  becomes 


3vA  _  .  3A  V3A3  ~  UpAR  , 

3x  Up  3x  VR  " 


Integrating  on  XO3/  and  noting  that  v(0,t)  =  0,  we  have 


v(x,t)  A( x, t ) 


-  u  A(x,t) 
P 


{V3A3  -  UpV 


V(x,t ) 
V 

R 


(23) 


or 


v(x,t)  A( x, t)  =  u  [Al  -  A( x, t ) ]  +  [v3A3  -  u  Ar]  -------  (24) 

P  P  R 

where  A^  =  AR  +  Ab¬ 
using  equation  (10)  in  (11)  we  obtain 


PLA 


3v 

3t~ 


f  PlAv 


3v 

3x 


(25) 


or,  integrating  on  x<s3 


P(x,t)  =  PQ(t) 


(x,t) 


dx'  . 


(26) 


The  objective  is  to  define  a  space-mean  pressure  using  the  pressure 
distribution  in  the  reservoir.  While  it  appears  feasible  to  simply  substitute 
equation  (24)  in  (26)  and  complete  the  required  integrations,  the  results 
would  be  complex.  Instead,  we  define  the  space-mean  pressure  on  [0,s-j]  as 

P(x,t)  dx  (27) 

where  lQ-j(t)  is  the  length  of  the  reservoir  from  0  to  s1  at  time  t  and  the 
velocity  v(x,t)  on  [0,3^]  is 


P(t)  = 


l01(t) 


: 


v  ( x#  t ) 


[v  A  -  u  A  "1 


X  ,  [O  i  S  1  ] 


Substituting,  we  obtain 


P(x,t)  *  PQ(t)  - 


T  °LX 


v3»3  - 


u  »„T 


”  v  ,A, 

-  U  A  1 

3 

.... 

"it” 

L  v°  J 

The  term  in  brackets  can  be  rewritten  as. 


V3A3  -  u^Ar  up  Al  (V3A3  -  up  Ar)  (V3  A3  -  up  Ar) 

Vi  ^  2  f  2 

R  VR 


V  ’  2 

R  v 


Using  equations  (29)  and  (30)  in  (27)  we  have 


,  o  V,A,  ”  u  A_ 

P0(t)  =R(t)  +-6“101(t>  - V~-~- 

R 


U  A  ) 


(u  +  v  )A  (v  A  -  u  A  ) 
£  3333  pR 

V 


IV.  EVALUATION  OF  INTEGRALS 


Consider  first  /  ^  v(x,t)  dx 


Equation  (24)  can  be  rewritten  in  the  form. 


v(x,t) 


Up  A(x) 


[v  A  -  u  A  ] 
_3  3  g  R 


_ L_ 

A(x) 


V(x) 


Integrating,  we  have 


dx 


+• 


(v  A  - 
3  3 


u  A  ] 

.je.5. 


V  (x) 
A(x) 


dx  • 


We  now  define  three  effective  lengths. 


L03(t> 


L03(t) 


1  f 8  3 

v  I  v(x) 

R  JO 

-ils-  f3  3  .Yiil. 

VR  Jo  A(X) 


dx 


L03  (t)  -  *L 


fa3  _dx 

Jo  *<«'»' 


and  let  iQ^ft)  represent  the  time-dependent  length  of  the  reservoir 
rear  of  the  liquid  chamber  to  the  exit  at  the  vent. 

Rewriting  the  integral,  we  have 


(32) 


(33) 


(34) 


(35) 


(36) 


from  the 


I 

I 


I  V(x,t 

J° 


)  dx  =  Up  [L-3(t)  -  l03<t)]  -  L^(t)  .  (37) 

U 


We  now  note  that 


fj3  *  ’  -X'fo 


v( x,t)  dx  +  u  v_ • 
P  3 


Using  equation  (37),  we  have 


s  a 

5  3  v(x,t)  dx  *  up  [t»Q3(t)  -  l03(t)  -  L^t)] 


A  [v  A  -  u  A  ] 

+  "3  T-  L03(t)  +  V3  *  — —X— «—  S>)  (39) 

L  L 


where  we  have  used  the  facts  that 


L03(t)  =  l01(t» 


-  u  ,  and  l.,(t) 
p  03 


-  u  ,  and  where  we  have  assumed  A  =  0 
P  3 


since  the  vent  area  is  variable  only  over  the  slant  section  of  the  bolt. 


Consider  now 


fs 

8t  Jo 


p  vA  dx  . 
Li 


Recalling  that  dp  /dx  =  0  ,  and  substituting  equation  (24),  we  have 


milt  WMMUJimwmMUPJM1 


UURIIIII 


i 


8. 


CL  I  3  *A(x.t)  ax  -  -fr  (PL[UpAL  a03(t)  -  ) 


♦  ‘v3R3  -  »PV  L03lt,1> 


(40) 


giving  us. 


i~  ■/: 


vA(x,t)  dx  = 


4P  <»A  a03(t) 


■>  -  A1.!1411 


*  "3  <RA  A  (V3A3  ‘  “pV  L03(t) 


A  fS 

*“  Jo 


vA  dx  - 


(41) 


V.  EVALUATION  OF  EFFECTIVE  LENGTHS 


In  order  to  simplify  the  effective  length  integrals  we  assume  that  the 
bolt  is  straight  and  consider  the  area  change  of  the  piston  only-  Then, 


LQ3(t) 


101  (t)  Al  1 

t  --  {112  [lQ1(t)  Al] 


2V 


it  1221222  1  4  4v ,  , 

4  ~r  ‘VW  «  3R1  -  rb  >  4  2  rb  (R1  'R2  1  -  12  (R1  -  R2  111 
M 


+  4;  {1H  [101(t)  \  +  V12]  +  [  S  A3  (113 


1 1 2 }  ]) 


(42) 


1  3 


L03  (t)  "  *  2  L01  (t)  ~v“ 


A  f  1  (t)  A 

i-  }+  -Ss-  {  -21--— 

R  VR  L  *M 


V(1/3)  R  2  -  r  2)  At  R,  +  r 

-1 . -2 . >  -Sr-  ^  <  t-  1  - 2  m  (  r~-;r  )] 

M  b  3  2  b 


rb  \ 

+  -----  m  (  —  ) 

3M  3 


R12  “  R22l  f  \  I"  123tl01(t)  \  + 

"  “’i?”*  J  L . a3 . 


+  ~2~  (  113  "  112  ’  +  101(t)  *23 


Lo3<t»  *  W*’  *  1  ln  <-r->  1  - 2  ln  1  -srrr-  *  1 

b  3  2  b 


23  A. 


L13(t)  =  L03(t)  ’  101(t) 


where  M  ■  (R1  -  R2 )  /  (x  -  ) 


While  the  effective  lengths  above  are  all  time  dependent,  the  time  derivatives 


do  not  contribute  terms  in  u  or  v,  • 

P  3 


»■» 


I  a 


VI.  PISTON  AND  LIQUID  ACCELERATION 


Substituting  (31)#  (39)  and  (41)  into  equations  (8)  and  (9)  we  have 


&p  (Mp  *  Pl\  L'o3M  -  »L*t,  (103(t)  - 


.Is. , . 

\  *  vR 


P,  (A  +  A  ) 
3  p  3 


•  ,  .  PLA3  *1  101<t> 

*  V3  1  PLA3  l03(t) . 6"»: - 


P  \  *  pLV3  A3 


i-  -l\  "o,11'  ♦ 


(u  +  v  )  A  (v  A  -  u  A  ) 

._J2 _ l _ l _ IA _ E_5_ 


+  pl(V3  A3  ‘  Up  V  L03(t)  + 


8pl  [‘ 

*  Jo 


vA ( x , t )  dx 


4P  [t03(t'  -  W*'  - 


.Is.  (t,  *  _I°ll!i.Is  1 

\  03  6vlt 


.  r  A3  ,2  A3  101(t) 

+  V3  1  ’a"  L03(t)  6V~  1 

L  R 


(P  -  P3) 


2  V3  "  hf  "  hf*  +  6  101(t) 


(u  +  v  )  A  (v  A 
£  3  3  3  3 

. 


-  u  V_  - 


v  A  -  u  A 
~3_3_„J2„R_  ,2 


L:.(t) 


Wu  WRWIWJWWW  V.VWWI 


Now  defining  G  and  F  as  the  right  hand  sides  of  (46)  and  (47),  we  have 

V  p  A  A  1^  (t) 

u„  [M  *  P-A  hi  At)  -  pA  (1__  (t)  -  ----  ) - - - ] 

p  p  L  R  03  L  L  03  A  6V 

L  R 

.  1  PLA3  AL  *01 *t3 

"  V3  tpLA3  L03(t) . 6V~~ - 3  =  G 

and 

3  AR  2  ^01  AR 

Up  [L03(t)  ~  103<t)  A  L03(t)  +  6V  3 

L  R 

A3  2  A3  X01(t) 

+  V3  [  *A~~  L03(t) . 6V~  3  -  F  ' 

L  R 

and  further  defining 

„  .  . 1  ,  plar  al  Loi(t)  , 

meff  “  1  pLA3  L03(t  6V~  3 

A3  2  A3  101^t3 

^ff  =  1  -A“  L03(t) . 6V~“  3 

A  V 

M  =  M  +  ----  m  _  -  PtAt  [1  ft)  -  ----  ] 

eff  p  A,  eff  KL  L  03  A 

3  Lt 

A 

L  =1  (t)  -  L3  (t)  +  — --  1 

eff  03  03  A3  eff 
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(48) 


(49) 


(50) 


(51) 


(52) 


(53) 


we  can  rewrite  equations  (48)  and  (49) 


u  M  .. 
p  eff 


v3  mef f 


v.  1 
3  eff 


u  L 
p  eff 


Solving  for  u  and  v_  ,  and  defining 
p  3 


M  *  -  M  -  p  A  [1  (t)  -  ] 

p  p  L  L  03  a. 


W**  -  L03(t) 


F*  «  - - F 

eff 


we  have 


•  ,  __1__  _“eff_ 

p  M  *  M  ■ 

P  P 


Jim.  .ah.  ,  .  g 

3  V  K"  V 


VII.  ADDITIONAL  EQUATIONS 


To  complete  the  description  of  the  liquid  chamber  we  need  equations  for 


the  piston  position  and  mass  flux. 


x  »  u 
P  P 


The  equation  of  state  for  liquid  monopropellants  is 


P  m  p  +  ---- 
R  °  K2 


1  3 


(62) 


where  the  density  is  determined  by 


m 

L 


(63) 


The  governing  equations  for  the  lumped  parameter  model,  equations  (58- 
63),  are  a  system  of  first  order  differential  equations,  which  have  been 
solved  in  standard  FORTAN  using  the  Adams  method  with  functional  iteration  in 
the  IMSL  library  DGEAR  on  an  IBM  PC-AT.  The  chamber  pressure  boundary 
condition  is  specified  as  an  input,  using  the  experimentally  measured  chamber 
pressure  from  a  2/3  charge  firing  of  a  Concept  VI  regenerative  liquid 
propellant  gun  test  fixture  at  the  Ballistic  Research  Laboratory. 


VIII.  RESULTS  AND  DISCUSSION 


Experimental  chamber  and  liquid  pressures,  and  piston  displacement  from 
the  2/3  charge  firing  used  in  specifying  the  chamber  pressure  boundary 
condition  are  presented  in  Fig.  2.  The  zero  in  time  has  been  chosen  as  the 
time  at  which  the  pressure  begins  to  rise  in  the  combustion  chamber  due  to  the 
influx  of  gas  from  the  igniter. 
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Figure  2.  Experimental  Chamber  Pressure,  Liquid  Pressure, 

Piston  Displacement 


The  piston  begins  to  move  at  about  1.25  ms,  travels  about  0.55  cm  and 


abruptly  stops  at  about  3.5  ms,  hestitates  briefly  and  then  again  accelerates 
and  smoothly  completes  its  stroke.  This  interrupted  piston  travel  is  a 
characteristic  of  the  Concept  VI  RLPG.  In  order  to  permit  the  piston  to  clear 
the  seal  on  the  nose  of  the  center  bolt,  the  piston,  liquid  reservoir,  and 
transducer  block  initially  move  rearward  approximately  0.55  cm  against  a  set 
of  Belleville  springs.  When  the  springs  are  fully  compressed,  the  transducer 


block  abruptly  stops,  as  does  the  reservoir  and  piston.  The  piston  then 
accelerates  rearward  again  as  liquid  injection  begins,  and  completes  its 


stroke. 


B 


1 


The  chamber  pressure  rises  steadily  to  about  20  MPa  in  about  2.4  ms,  and 
then  maintains  that  value  until  about  4  ms  when  liquid  propellant  combustion 
begins.  The  chamber  pressure  then  rises  smoothly  to  its  maximum  value,  drops 
slightly  as  the  piston  reaches  the  rear  taper  of  the  bolt  at  about  8  ms,  and 
then  drops  sharply  at  burnout  as  the  piston  completes  its  stoke  at  about  9  ms. 

The  propellant  in  the  reservoir  is  much  stiffer  than  the  combustion 
gases,  and  thus  reflects  the  abrupt  variations  in  piston  motion.  As  the 
Belleville  springs  begin  to  compress,  a  small  oscillation  in  liquid  pressure 
is  observed  at  about  3  ms.  When  the  transducer  block  suddenly  stops  at  about 
3.5  ms,  the  momentum  of  the  piston  is  absorbed  by  the  liquid,  producing  the 
relatively  large  hydraulic  pressure  oscillations  from  3*5  ms  to  5.0  ms. 
Initially  these  oscillations  are  undamped;  however,  as  the  injection  area 
opens,  the  oscillations  are  rapidly  damped.  Similarly,  as  the  piston  reaches 
the  rear  taper,  which  reduces  the  liquid  injection  area,  the  liquid  pressure 
rises  sharply  as  the  piston  is  decelerated.  The  liquid  reservoir  gage  fails 
just  as  damping  begins  at  about  8  ms. 

No  attempt  has  been  made  here  to  simulate  the  motion  of  the  transducer 
block  against  Belleville  springs.  Instead,  the  zero  in  time  is  chosen  to  be 
the  point  at  which  the  piston  stops  due  to  completion  of  transducer  block 
motion,  and  the  :nitial  conditions  for  the  solution  of  the  ordinary 
differential  equations  are  then  taken  from  experimental  data. 


A  simulation  was  first  made  using  Equations  (1-3)  as  a  baseline.  A 
description  of  the  input  data  can  be  found  in  Appendix  A.  The  results  show 
good  agreement  in  the  predicted  and  measured  liquid  pressures.  A  comparison 
of  the  input  discharge  coefficient.  Equation  (3),  the  calculated  discharge 
coefficient.  Equation  (64),  and  the  experimental  discharge  coefficient 
determined  by  Coffee  is  presented  in  Figure  3. 
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Figure  3.  Discharge  Coefficients  from  a  Simplified  Flow  Model 

Without  Inertial  Terms 


The  steady  state  discharge  coefficient  is  calculated  from  Equation  (3). 
It  varies  slightly  due  to  the  variation  in  Reynolds  Number  over  the  ballistic 
cycle,  but  this  variation  is  quite  small. 


The  predicted  discharge  coefficient  initially  exhibits  large 
oscillations,  which  are  rapidly  damped.  These  oscillations  are  a  direct 
reflection  of  the  hydraulic  oscillations  in  the  propellant,  and  the  fact  that 
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the  liquid  flow  in  the  orifice  cannot  adjust  to  the  rapid  fluctuations  in 
liquid  pressure.  The  result  is  an  artificially  large  (and  totally  ficticious) 
"discharge  coefficient"  calculated  from  Equation  (64).  After  the  oscillations 
damp  out  at  about  1.5  ms,  the  predicted  discharge  coefficient  rises  slowly, 
reaching  the  steady  state  value  at  about  3.5  ms  as  the  pressures  reach  the 
steady  state  operating  regime  near  maximum  pressure.  The  "rise"  in  the 
predicted  discharge  coefficient  coincides  with  the  rapid  rise  in  chamber  and 
liquid  pressures,  and  the  acceleration  of  the  piston  to  its  steady  state 
velocity.  This  indicates  that  even  in  the  simple  model,  the  liquid  velocity 
lags  the  pressure  drop  from  the  reservoir  to  the  chamber,  resulting  in  an 
apparent  "low"  value  for  the  discharge  coefficient. 


The  experimentally  determined  discharge  coefficient  is  quite  different  in 
magnitude  from  the  predicted,  but  displays  some  general  similarities.  The 
early  oscillations  are  present,  though  of  a  different  frequency,  and  reduced 
both  in  magnitude  and  mean  value.  The  increase  in  the  experimental  discharge 
coefficient  to  its  maximum  value  occurs  over  the  same  interval  as  in  the 
predicted  case.  However,  the  experimental  discharge  coefficient  begins  at  a 
much  lower  value  than  predicted,  about  .04  ms,  and  peaks  at  about  1.1.  The 
mean  value  of  the  experimental  discharge  coefficient  over  the  steady  state 
interval,  3.0  to  5.0  ms,  is  about  0.95,  in  comparison  to  a  predicted  value  of 
about  0.85.  Thus  the  simple  lumped  parameter  simulation  is  not  only  unable  to 
reproduce  the  transient  behavior  of  the  discharge  coefficient,  but  is  also 
unable  to  account  for  the  high  mean  values  of  the  experimental  discharge 
coefficient. 


The  effects  on  the  predicted  discharge  coefficient  of  extending  the 
control  volume  to  include  the  entire  propellant  reservoir  are  shown  in  Fig. 
4.  The  experimental  discharge  coefficient,  the  predicted  discharge 

coefficient  using  Equations  (1)  and  (2),  and  the  predicted  discharge 
coefficient  from  the  model  developed  here  are  presented.  The  two  predicted 
discharge  coefficients  are  very  similar  in  structure.  Both  show  the  large 
oscillations  discussed  above,  and  both  display  the  slow  increase  as  the 
systems  approaches  steady  state  operation.  However,  the  magnitude  of  the 
discharge  coefficient  obtained  using  our  model  is  significantly  higher  than 
that  for  the  simpler  model.  Our  predicted  discharge  coefficient  agrees  quite 


well  with  the  mean  value  of  the  experimentally  derived  data  in  the  region  of 
steady  state  operation.  The  lack  of  agreement  in  the  early  values  of  the 


discharge  coefficient/  and  the  rise  to  steady  state  is  not  apparent.  However * 
it  appears  that  there  is  some  uncertainty  in  the  actual  piston  position,  which 
could  significantly  affect  the  computed  injection  area,  and  thus  the 
experimental  discharge  coefficient.  We  discuss  the  uncertainity  in  the 
initial  piston  position  in  greater  detail  below. 
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i 

|  Figure  4.  Discharge  Coefficients  from  Model,  Gough  Formulation, 

I  and  Experiment 


'ttie  predicted  and  experimental  liquid  pressures  and  piston  displacements 
are  presented  in  Figs.  5  and  6  respectively.  The  predicted  and  experimental 
liquid  pressures  agree  quite  well  over  most  of  the  ballistic  cycle,  however, 
the  predicted  pressure  exceeds  the  experimental  value  by  about  10%  at  peak 
pressure.  There  is  an  obvious  discrepancy  in  the  predicted  piston  motion. 
The  experimental  data  indicate  that  the  piston  comes  to  a  stop  when  the 
transducer  block  completes  its  motion,  and  briefly  remains  motionless  before 
again  accelerating.  In  comparison,  the  predicted  piston  motion  is  only 
slightly  perturbed,  and  the  piston  continues  to  move  with  little  hesitation. 

In  the  case  of  the  initial  piston  motion,  the  predicted,  hydraulic 
pressure  oscillations  are  damped  more  rapidly  than  in  the  experimental  data. 
There  is  some  uncertainty  in  the  actual  displacement  of  the  transducer  block 
and  thus  the  piston.  The  piston  is  very  close  to  the  position  where  the 
injection  area  opens  when  the  transducer  block  completes  its  stroke. 
Therefore,  a  slight  discrepancy  in  the  piston  position  can  have  a  large  effect 
on  the  injection  area  and  thus  system  damping.  The  comparison  of  simulated 
and  experimental  data  would  suggest  that  the  initial  piston  position  assumed 
in  the  simulation  is  incorrect,  and  that  the  injector  is  opened  too  early  in 
the  process.  As  a  result,  the  piston  completes  its  stroke  earlier  in  the 
simulation  than  in  the  actual  experiment. 

In  the  case  of  the  maximum  liquid  pressure,  the  discrepancy  is  related  to 
the  problem  with  the  piston  motion.  In  the  experimental  data,  the  piston 
reaches  the  rear  taper  after  P^ax'  anc*  the  chamber  pressure  shows  an  almost 
immediate  drop  as  the  injection  area  begins  to  decrease.  However,  in  the 
simulation,  the  piston  reaches  the  rear  taper  near  PMAx'  as  experimental 
chamber  pressure,  which  is  an  input,  is  still  increasing.  The  combination  of 
a  decreasing  injection  area  and  an  increasing  chamber  pressure  would  be  an 
increase  in  liquid  pressure. 


IX.  CONCLUSION 


A  lumped  parameter  model  for  the  motion  of  the  regenerative  piston  and 
the  injection  of  liquid  propellant  has  been  presented.  This  model  includes 
the  entire  propellant  reservoir,  and  thus  the  effect  of  the  inertia  of  the 
propellant  in  the  reservoir  on  the  process.  The  coupling  between  piston 
motion  and  liquid  injection  are  fully  included  as  well. 

The  results  of  computer  simulations  using  this  model  are  compared  both 
with  a  somewhat  simpler  model  and  with  data  derived  from  experimental  gun 
firings.  The  results  of  these  comparisons  are: 


1.  The  simpler  lumped  parameter  model  exhibits  a  transient  behavior 
similar  to  that  of  the  experimental  data,  but  the  magnitude  of  the  discharge 
coefficient  is  incorrect. 

2.  The  early  oscillations  are  related  to  hydraulically  induced  pressure 
oscillations  in  the  propellant  reservoir,  and  do  not  represent  a  real 
variation  in  the  discharge  coefficient. 

3.  The  slow  rise  in  the  discharge  coefficient  from  3.0  ms  to  5.0  ms 
corresponds  to  the  period  of  rapid  piston  acceleration  and  rapid  pressure 
rise.  During  this  period,  the  discharge  coefficient  remains  somewhat  below 
the  steady  state  value.  This  appears  to  be  due  to  the  injection  velocity 
lagging  the  pressure  drop  during  the  rapid  approach  to  steady  state. 

4.  The  maximum  value  of  the  magnitude  of  the  discharge  coefficient 
obtained  from  the  simpler  model  is  significantly  less  than  the  mean  value  of 
the  experimental  data.  The  cause  of  this  discrepancy  is  not  apparent. 

5.  The  discharge  coefficient  obtained  from  a  simulation  using  the  model 
developed  here  exhibits  transient  behavior  almost  identical  to  that  for  the 
simpler  model.  However,  the  magnitude  of  discharge  coefficient  is 
significantly  higher  than  in  the  case  of  the  simpler  model,  and  agrees  quite 
well  with  the  mean  value  of  the  experimental  data  in  the  steady  state  region. 

6.  There  are  discrepancies  between  the  experimental  liquid  pressure  and 
piston  motion  and  the  simulation  using  the  model  presented  here.  It  appears 
that  these  discrepancies  are  the  result  of  uncertainties  in  the  initial  piston 
position. 

The  next  objective  in  this  study  is  correction  of  the  discrepancies. 
Future  work  will  focus  on  the  elimination  of  simplifying  assumptions  involving 
the  system  geometry,  and  the  investigation  of  a  more  advanced  RLPG 
configuration  in  which  the  Belleville  springs  are  eliminated.  The  Lagrange 
pressure  distribution  will  be  extended  to  include  the  injection  orifice,  and 
the  resulting  pressure  and  velocity  distributions  will  be  compared  with  the 
results  of  one  and  two  dimensional  simulations. 
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Cross-sectional  area  of  the  vent,  cm^ 

Cross-sectional  area  of  the  liquid,  cm^ 

2 

Cross-sectional  are  of  the  piston  on  chamber  side,  cm 
Cross-sectional  are  of  the  piston  on  reservoir  side,  cm^ 

Discharge  coefficient  for  the  liquid  into  the  chamber 
Diameter  of  the  hole,  cm 

Friction  loss  of  the  liquid  in  the  vent 
Entrance  loss  of  the  liquid  into  the  vent 
Bulk  modulus  at  zero  pressure,  MPa 
Derivative  of  bulk  modulus,  MPa 
Length  of  the  liquid  column  from  sq  to  s^,  cm 

Length  of  the  liquid  column  from  sQ  to  s^,  cm 

Length  of  the  liquid  column  from  s,)  to  s2,  cm 

Length  of  the  liquid  column  from  s1  to  s^,  cm 

Length  of  the  vent 

Mass  of  the  piston 

Mass  flux  of  liquid  into  the  combustion  chamber,  gm/s 

Breech  pressure ,  MPa 

Combustion  chamber  pressure,  MPa 

Liquid  reservoir  pressure,  MPa 

Space-mean  pressure  in  the  reservoir,  MPa 

Radius  of  the  bolt  at  x^,  cm 

Radius  of  the  bolt  at  x2,  cm 

Radius  of  the  bolt  at  x,  cm 

Radius  of  the  piston  at  s^,  cm 

Radius  of  the  piston  at  a^,  cm 
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Volume  of  liquid  from  8^  to  s^,  cm"* 
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Volume  of  the  liquid  reservoir,  cm 
Position  of  the  piston,  era 
Density  of  the  liquid  initially,  gm/cm 
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TABLE  A-1.  Input  Data  for  Simulation 


RLF'LCH — ROUND  3 — 30MM 


COMBUSTION  CHAMBER  AREA  -  44.84700 
PISTON  AREA — C  CH  SIDE  -  34.32600 
PISTON  AREA — RES  SIDE  *  23.27800 
LENGTH  L  PRIME  «  1.43200 
LENGTH  OF  VENT  «  1.04000 
PISTON  MASS  »  2109.20000 
VOLUME  LIQUID  -  172.63196 


VENT  OPTION 

STRAIGHT  LENGTH  OF  PIST 
MAX  PISTON  TRAVEL 


2 

5.94680 

7.37S80 


DENSITY  LIQUID 

Kl 

K2 

INLET  LOSS 


1 . 43700 
5661. 10000 
9.26490 

0.62000 


FRICTION  LOSS  OPTION 
FRICTION  LOSS 


0.00000 


TIME-C  CH  PRES  DATA  FILE:  A: PT0FF64. DAT 


GEOMETRY  DATA  FILE: 


A: QFF544Z . DAT 


GRAPH  DATA  FILE: 


A i 0F544IM.GR A 


INITIAL  PR  IN  RESERVOIR  =  29.00000" 

INITIAL  VEL  IN  VENT  -  0.00000 

INITIAL  PISTON  VELOCITY  -  358.00000 

INITIAL  PISTON  POSITION  =*  0.00000 

INTEGRATOR — TINC  «  0.00010 

INTEGRATOR — EPS  -  0.00001 

INTEGRATOR — METH  -  1 

INTEGRATOR — MITER  »  0 

INTEGRATOR — KWRITE  »  O 


DIFFERENTIAL  EQUATION  SET:  1 
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